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FEEDBACK CONTROL FOR DAMPING 
A SYSTEM OF LINEAR OSCILLATORS 

ALEXANDER OVSEEVICH AND ALEKSEY FEDOROV 


Abstract. The paper is devoted to a design of a common bounded feedback 
control steering a system of an arbitrary number of linear oscillators to the 
equilibrium. At high energies, the control is based on the asymptotic theory 
of reachable sets of linear systems. With decreasing of the energy, a similar 
control with a reduced upper bound is used. On the final stage, the control 
is constructed by using the method of common Lyapunov functions. Special 
attention is paid to the cases of one and two oscillators. 
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1. Introduction 

In control theory, a class of problems related to the steering of a system from an 
arbitrary initial state to a given manifold in minimum time is well studied. Strictly 
speaking, this problem consists of two. The first of them is to bring the system to 
a given manifold. The second problem requires minimization of the motion time. 
There are two main approaches to design of the corresponding control. First, a 
feedback control (synthesis), where the control u{x) is a function of the current 
state of the system. Second, a feedforward control u{t), which is a function of the 
elapsed time. One of the classical results of control theory is the analytic design of 
the feedback control for a single linear oscillator based on the Pontryagin maximum 
principle [1]. 

1.1. Problem statement. In the present paper, a next in complexity problem is 
considered: 


Problem 1. Design a feedback control for steering a system of N linear oscillators 
with eigenfrequencies uJi 

(1.1) X = Ax ^ Bu^ X = (xi, ^ 1 ,..., xat, ^at)* G rt G M, \u\ < 


( 1 . 2 ) 


0 1 

—uol 0 


to the equilibrium in minimum time. 




1 
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/ 0 \ 
1 


0 

V1 / 


Here and in what follows * stands for the transposition. 

It looks like that an analytic construction of the optimal feedback control by 
using maximum principle methods in Problem 1 is impossible, and even the search 
for a numerical solution is a hard problem. 

A mechanical model for system (1.1)-(1.2) is a system of N pendulums with 
eigenfrequencies uJi attached to a cart moving with bounded acceleration where 
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a) 





b) 


Figure 1. A mechanical model of the system: a) system of N 
pendulums with eigenfrequencies uji attached to a cart; b) masses, 
attached via springs to a body, which moves with acceleration. 


the vertical deviations Xi of pendulums are components of the state vector (Fig. 
la). In the other interpretation, components of the vector are the displacements Xi 
of masses, attached to a body moving with acceleration u (Fig. lb). 

1 . 2 . Proposed approach. A non-optimal feedback control, which steers the sys¬ 
tem to the equilibrium, will be used. The control is asymptotically optimal: a ratio 
the time motion to zero by using this control to the minimum one is close to one if 
the initial energy of the system 

N 

(1.3) E = 1/2^ {xi‘^ + 

i=l 

is sufficiently large. Proof of the asymptotic optimality of the suggested control is 
beyond our paper and it will be presented elsewhere. 

In our paper, an assumption of absence of resonance, ie., non-trivial relations 
between eigenfrequencies of the form, 

N 

(1.4) rriiUJi = 0, where 0 7 ^ m = (mi, • • •, ^at) G 

i=l 

is used. In particular, this implies that for system (1.1)-(1.2) the Kalman control¬ 
lability condition holds [2]. The suggested control works also in the resonant case, 
however its quasi-optimal properties are lost. 

To solve the problem, three strategies are used in sequence. 

1 . High energy zone. At high energies, the normal vector to the approximate 
reachable set, close to the exact one in the long run, serves as a momentum 
p{x) (Fig. 2a) [3-5]. 

2 . Intermediate zone. The obtained control is applicable also at low energies, 
but its quasi-optimal properties in this case are lost. First, this control af¬ 
fects the system like a dry friction, so that in some states it prohibits any 
motion. Second, the motion may occur in a small neighborhood of a limit set 
(attractor), which does not contain the equilibrium position. Application of 
the control with reduced upper bound (second stage of the control) allows 
to postpone this undesirable pulling into the attractor. 

3. Low energy zone. At the final third stage, an approach to the design of a 
local control in the feedback form based on the common Lyapunov functions 
is used. This approach works in a small neighborhood of the equilibrium 
position (Fig. 2b). To reach this small neighborhood, it should necessarily 
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Figure 2. Structure of the suggested control: the region bounded 
by the solid line is the region in which the control based on the 
asymptotical theory of reachable sets is used; the region bounded 
by the dashed line is the region in which analogical control with 
reduced upper bound is used; the region bounded by the dotted 
line is the region in which the local feedback is applied. The point 
O is the equilibrium position. 



contain the basin of attractor of the preceding control. The reduction of the 
basin achieved at the previous stage makes this objective possible (Fig. 2c). 

It should be noted that there exist other methods of control design in the feedback 
form, e.^., based on the Kalman approach to the feedforward control for linear 
systems [9,10]. The cited works also include estimates for the motion time under 
the control, which are comparable with the optimal time in the sense that the ratio 
the motion time under the the proposed control to the minimum one is bounded. 

1.3. Goal of the work. A possibility of a design containing three steps for the 
system of linear oscillators has been demonstrated for the first time in Ref. [11]. 
A brief summary of results on the asymptotic optimality of the suggested control, 
partially based on a new lemma on the observable time-invariant linear systems, 
and on the study of the character of singular motion, can be found in Ref. [12]. 

In the present paper, we emphasize the computational and “algorithmic” aspects 
of the suggested control. The paper is organized as follows. Section 2 is devoted 
to a design, computation, and properties of the control used in high energy and 
intermediate zones. We give the explicit expressions for the support function of 
the reachable set of system (1.1)-(1.2). In the general case, the support function is 
expressed via the Gelfand hypergeometric function. In the case of two oscillators it 
is expressed via the elliptic integrals. Methods for solution of the master equation 
determining the control are indicated. In the case of two oscillators, the method 
is based on the solution of a transcendental equation in elliptic functions. In the 
case of an arbitrary N^ we state an optimization problem such that its solution is 
equivalent to the solution of the master equation. For the case of two oscillators 
N = 2 the solution of the optimization problem is equivalent to the solution of a 
transcendent equation in elliptic functions. Section 3 is devoted to construction of 
the control in a region close to the terminal point. In particular, common Lyapunov 
functions for A" = 2 and N = 1 are calculated, it is shown also that for the arbitrary 
N the common Lyapunov function is an even integer-valued matrix. In section 4, 
main matching relations for the control at different stages are presented. In Section 
5, the suggested strategy is illustrated for the case of a single linear oscillator. 

Proofs of all the statements can be found in [13]. 
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2. Control at high energies 


2 . 1 . Asymptotic theory of reachable sets. Being simultaneously an instrument 
and an object of studies, reachable sets occupy a central place in control theory. 
Therefore, their study is a fundamental problem [14]. Recall that the reachable set 
V{T) is defined as the set of points reachable from zero at time T. 

One of the main result of Ref. [4] as applied to the considered system ( 1 . 1 )-( 1 . 2 ) 
is as follows: the reachable set V{T) has as T ^ oc asymptotic form Til, where 
is a fixed convex body. The result is stated more precisely in the following theorem. 


Theorem 1. Suppose that the momentum p is written in the form p = 
where is the dual variable for Xi, rji is the dual variable for yi, and Zi = {yf + 
In the ease of absenee of resonanee (nontrivial relations between eigen- 
frequeneies of the form (l.f)), the support funetion Ht of the reaehable set V(T) 
has as T ^ oo the asymptotie form 
( 2 . 1 ) 


Ht{p) = 


T 

(2^ 


n27T ^277 

Jo Jo 


N 

E 


Zi cos ifi 


dipi... d(fiN + o(T) = TS^{z) + o(T), 


and the support function of the compact set Q, is given by the main term Sj(z). 


The support function of any subset M € is defined as 


( 2 . 2 ) 


Hm(0 = sup(^,a;), 

xeM 


where angle brackets stand for the scalar multiplication in 


2 . 2 . Approach to the control design. Geometrically, the Pontryagin maximum 
principle says that the momentum (vector of adjoint variables) t/; in a point x is the 
inner normal to the reachable set V{T) {T(x) being the time of reaching x from 
zero). 

The idea of our approach to control design is to use the set TVt as an approxi¬ 
mation to the reachable set V{T), and the normals to TO as momenta. If the phase 
vector X G V = lies at the boundary of the set TO, then 


(2.3) 


T-x= 


for a certain momentum p = p{x). We note that the support function is differ¬ 
entiable at p 7 ^ 0 and Eq. (2.3) has a unique solution because of the smoothness of 
the boundary of the set O proved in Ref. [5]. The feedback control is given by 


(2.4) 


u{x) = —sign(T,p(x)). 


In what follows, we will also use the control in the form 


(2.5) 


uu{x) = Uu{x), \U\ < 1. 


2.3. Computation of the control. In the coordinates notation, Eq. (2.3) has the 
form 

(2.6) T-\xi,yi) = zp , i = l,...,N, 

where Zi = ( 77 ? + In order to solve Eq. (2.6) we calculate first the point 

3 of a sphere with positive-homogeneous coordinates {zi : • • • : zn)- Here the 

sphere is regarded as a set of directions of non-zero vectors in To this 
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end, define an “energetic” vector e = (e^) G where and 

obtain from Eq. (2.6) that 

(2.7) T-^ei = j^, i = 

dzi 

The solution of Eq. (2.6) defines an inversion of a map from one 2A/'-dimensional 
manifold to another one, while the solution of (2.7) reduces to inversion of map 
of (A^ — l)-dimensional manifolds. Like Eq. (2.6), the Eq. (2.7) possesses a unique 
solution [5], which can found as follows. 

In the case N = the support function has the form S){z) = ^\z\, and the 
solution is 3 = e, while the control has the dry-friction form u = —signal. If N = 2, 
then Eq. (2.7) is one-dimensional: the function S) can be expressed through elliptic 
integrals, point (ei : 62 ) G defines problem (2.7) completely, and the solution 
3 G 5^ can be also found by solving a transcendental equation in elliptic functions. 


2.3.1. Formula for the support function. The support function of the convex body 
Cl is given by the main term of the asymptotic expression in ( 2 . 1 ) 

N 

.Zi cos (fi 


( 2 . 8 ) 


Hq{p) =^{z) = J 


i=l 
2 I 


dif^ where 2: = (zi,..., zn) G 




Eor = 1, we obtain S){z) = in the case N = 2 the function 
(2.9) ^{z) = J \zi cos (fi -h Z 2 COS (P2 I dip 

can be expressed via elliptic integrals (we will dwell on this subject below). In the 
general case. 


( 2 . 10 ) 


Sj{z) 


(27r)^ / 




N 






N 




i=l 


is an Euler-type integral that defines a generalized hypergeometric function in the 
sense of I.M. Gelfand 


2.3.2. Case of two oscillators. Our previous arguments imply that the case N = 2 
is special in some sense, because in this case one can solve Eq. (2.3) by a simpler 
method. Indeed, in this case 

277 

COS 0isign(2(i cos (pi -h Z2 cos (p2)dip. 

Eor instance, consider the index i = 1 and make the inner integration wrt ip 2 . 
Taking into account the non-negativity of Z 2 , one needs to compute the integral 

1 2 

( 2 . 11 ) — / sign(—C + cos(/:? 2)<^^2 = — arccosC — 1 , |C| < 1 , 

Z7T Jo TT 

where C = /ccos^i, k = —zilz 2 . It is obvious that by permutation of the indexes, 
we can assume that \k\ < 1. Note that the latter assumption breaks the symmetry 
between zi and Z 2 . At \k\ < 1 , from ( 2 . 11 ) we obtain that 


^ ^ ^ 7 

dzi 271 Jo 


( 2 . 12 ) 


dS) 

dzi 


1 7^^ 

Jo 


-— = ^ / cosipi diTCCOs{kcosipi)dipi, 


because cos (pidcpi = 0. After integration by parts, integral ( 2 . 12 ) can be written 
in the elliptic form 


(2.13) 


n2TT 

L ' 


cos if arccos(7: cos ip)ip 


n27r 

Jo 


k sin^ (p 




COS^ (p 


dip 
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and finally we get for the following formula for the derivative of the support function 

k sin^ cp 




y/l — k‘^ COS^ ip 


dip, where k = —zxjz^^ 


which holds true provided that \k\ < 1. For computation of we need to compute 
the inner integral 

1 2 

— / cos (/:? 2 sign(—C cos ip 2 )dip 2 = — sinarccosC, if \C\ < 1 , 

Z7T Jq TT 

wherefrom we obtain 

(2J5) 


Jo 


dz2 

Note that the asymmetry between integral formulas (2.14) and (2.15) is illusive: 
the change of variables zi ^ 2:2 implies the change of parameters k ^ k~^. Under 
this change, the integrals 

p27T 1 • 2 p27T 

(2.16) Ii{k) = / sm ip l 2 {k) = / ^/Y^^k?co^^dip, 

Jo y 1 — k‘^ cos‘^ ip Jo 

regarded as meromorphic functions of k, are transposed: Ii{k~^) = l 2 {k). 

The key equation (2.7) that defines control (2.4), has the form of equation for 
k = -Z 1 IZ 2 : 


(2.17) 


? = 
ei h 


For numerical experiments, it is useful to employ a representation for (2.12) and 
(2.15) via canonical Legendre elliptic integrals of the first and second kinds 

(2.18) Aik) = pl-k'^S ’ ^(y) = j yi simpdip, 

where E {y) is the complete Legendre integral of the second kind, 

(2.19) 


Bik) = - 


yjl — kP' 


k^ 


k‘^ 


-K 




k‘^ 


>C{y)=r 

J 0 


/2 


dip 


— y‘^ sinip^ 


where JC {y) is the complete Legendre integral of the first kind. 
Then expression (2.17) can be rewritten in the form 


F(k) = ^ 

^ ’ dz^j dzi Aik)' 

A 

(4 - 4)dF 


( 2 . 20 ) 


The graph of the function F{k) is presented in Fig. 2 . Note that the support function 
has the form 


1 ( 7 ^ - 7 ' 

( 2 . 21 ) ^i,,^,,)= -7== 

^ Jo V ^2 “ 


^ COS^ ip 


kl| < k 2 |- 


Note that the function dS)ldzi can be regarded as a function 5 * = 5 ^( 3 ) on 
Assuming that gi are known, we obtain the final formula for the momentum 

(2.22) i^i,Vi) = ^i^iXi,yi), i = l,...,N. 

If we know the point 3 = ( 2)1 : • • • : zn) ^ then due to the non-negativity of 

the coordinates Zi the direction of the momentum p{x) is defined by formula ( 2 . 22 ) 
uniquely: the unknown positive pre-factor T~^ does not play a role. Control (2.4) 
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F{k) 



k 


Figure 3. Plot of the function F{k) governed by Eq. (2.20). 


depends on the direction of the momentum only, and therefore it can be efficiently 
found. It has the form 


(2.23) 


N 


u{x) = -sign y] Qi 


-1 


ZiVi 


\i=l 


2.3.3. Case of an arbitrary number of oseillators. For an arbitrary number of linear 
oscillators, thanks to the Kuhn-Tucker theorem the search for a solution of Eq. (2.7) 
is equivalent to solving the following optimization problem 


(2.24) (6:^) ^ max, provided that S){z) < 1. 

It is obvious that the restriction S){z) < 1 is equivalent to the restriction S){z) = 1. 
The hypersurface {S){z) = 1} is strictly convex due to the relation 


(2.25) 




VjCOS(^i da{(p), 


where the integration is over V = {z E : f{z) = O} and 


(2-26) f{z) = {2TT)^^df'^^ 

is the canonical volume element on V. 

It follows from relation (2.25) that if the vectors f and z are not collinear, then 
(2-27) (0{z)C,f^>O. 

However if the vector ^ is tangent to the hypersurface {S){z) = 1} in the point 
then these two vectors cannot be collinear. Indeed, the tangency condition 
{dS)ldz^ z) = 0 contradicts the Euler relation {dS^jdz^ z) = Sj(z) >0 for the homo¬ 
geneous function Sj. Strict convexity of the hypersurface {S){z) = 1} means unique¬ 
ness of the solution of the optimization problem (2.24). From this strict convexity 
it follows that the function f = is strictly convex as well. 

In the same time, the optimization problem (2.24) is equivalent to 


(2.28) 


(e^z) max, provided that f{z) < 1. 
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if zi 7 ^ Z 2 are solutions of (2.28), then (e, zi) = {e, Z 2 ) and f{zi) = 1. This means, 
however, that 

(2.29) (e, 54 ^) = (».-'i) aiid / < 1, 

which contradicts the optimality of vectors Zi. 

Thus, the equivalent optimization problems (2.24) and (2.28) can be solved by 
using the well-known efficient methods, e.g., using the Matlab Optimization Tool¬ 
box. 


2.4. Asymptotic optimality of the control. Define a polar-like coordinate sys¬ 
tem, which is well suited for description of the motion under the control u. Write 
the phase vector x in the form x = where p > 0 and uj = dCl. In the terms of 
Eq. (2.3), p = T, and 0 = dHo^/dp. In these coordinates, the equation of motion 
has the form 


(2.30) p 



0 = + 


1 

P 


^Bu + (j) 



If A = 1, we get the proper polar coordinate system in the plane. 

For the function p = p(x), an eikonal-type equation holds 

(2.31) iJn(p) = l, 

It is dual to the equation p{dHQ/dp) = 1 of the surface uj. Eq. (2.31) can be used 
for averaging the right-hand side of the first identity in (2.30) with respect to time, 
and the proof of the following statement about the asymptotic optimality of control 
(2.4) is based on the equation: 


Theorem 2. Consider the evolution of the value p under eontrol (2.4) • Let M = 
min{p(0), p(T), T}. Then as M ^ +oc we have 

(2.32) (p(0)-p(T))/T = 1+0(1). 

Under any other admissible eontrol, 

(2.33) {p{0)-p{T))/T<l+o{l). 

The proof of this statement is given in Ref. [13]. 


2.5. Comparison with the maximum principle. Solution of the linear minimum¬ 
time problem reduces completely to the boundary value problem for the Pontryagin 
maximum principle corresponding to the Hamiltonian 


(2.34) h{x,fj) = {Ax,fj) -h — 1 = max{(Ax,'0) -h {Bu,fj) — 1}, 

where | • | is the Euclidean norm, and the maximum is taken over the interval 
{r G M : |r| < 1}. The problem has the form 


(2.35) 


X = Ax -h Bu, fj = —A'^fj, 

u = sign(H,'0), x(0) = xo, x{T) = 0, h{x,fj) = 0. 


One can approach the issue of optimality of control (2.4) by comparing the dif¬ 
ferential equations of the motion under the control with equations of the Pontryagin 
maximum principle. To this end, one need to understand how the momentum p{x) 
in the Eq. (2.4) changes with time. This description is given by the following equa¬ 
tion 

(2.36) p = —A*p + Bu, where B = 

Note that if the second term Bu is absent from the latter equation, then the equa¬ 
tion for fj = —p would coincide with the equation of the maximum principle for 
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o2 

the conjugate variables. However, the matrix ^ is a homogeneous function of x 
of degree — 1 , and, then the mentioned second term has the order of magnitude 
0{l/\x\) at large x so that it is small. Also note that the condition of maximum 
u = sign(H,'0) = —sign(H,j 9 ) holds for control (2.4). The motion under control 
(2.4) is governed by the Hamiltonian H = {Ax, ip) + \ {B,'ip) \ — \{B,dp/dx)\, which 
is in some sense close to the Pontryagin Hamiltonian h{x,ip). The difference be¬ 
tween the Hamiltonians 1-L and h is 1 — | {B, dp/dx)\. The average value | {B, dpjdx)\ 
is close to one for sufficiently large x, which follows from Theorem 2. 

Thus, for the vector (t,' 0), where = —dpjdx the equation of the maximum 
principle holds “on average” with a small error at large x. 

3. Feedback near the terminal point 

The idea of a design of local feedback control, which is a base of the third stage 
of the suggested control, goes back to [ 6 ] and uses a preliminary reduction of system 
(!.!)-( 1 . 2 ) to a canonical form by using transformations 

(3.1) A i-A A + HC, u^u — Cx, A^D~^AD, B^D~^B, 

that correspond to adding a linear feedback control, and to coordinate changes. 
Formulate the result as a Lemma. 

Lemma 1. By using transformations (3.1), system (1.1)-(1.2) reduces to the form: 

(3.2) ? = % + ^u. 


(3.3) 

21 = 

/ 0 
-1 0 
-2 0 


, ® = 

0 

0 



V 

-2N + 1 0 J 


\oJ 


Wherein the matrix of the linear feedback has the form 
(3.4) C = (ci OC 2 O ... CivO), Ck = 


ijtk 

The gauge matrix D has the following form. Define 2x2 matrices 

1 \ 

, where Xk = ^ 


(3.5) 


d,, = (-i)^-'Ar' 


(2(i-l))! 


0 




i^k 


Then, D is the N x N matrix {dij) of 2x2 blocks dij. 


As an existence theorem of a canonical form without explicit formulas for ma¬ 
trices C and D, Lemma 1 is a particular case of the Brunovsky theorem [15]. By 
following [ 6 ], we introduce a matrix function of time related to system (3.3): 

(3.6) 5(T) = diag(T\T2^...,T^^)-^ 

In what follows the parameter T will become a function T = !!(}:) of the phase 
vector. Define the matrices in accordance with [6,8] 

q = (qij), qij = p - x)dx = [{i + j){i + j - 1)]"\ 

(3.7) 

0 = q-i, £=-1^*0, fm = diag(l,2,...,2Af). 

Define the feedback control by the equation 

(3.8) u(y) = mnWh 
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where the function 1 = T(j:) is defined implicitly by the following relation: 

(3.9) = 

A basic result on the steering of the canonical system (3.2)-(3.3) to zero is as 
follows: 


Theorem 3. Consider system (1.1)-(1.2) in the eanonieal form (3.2)-(3.3) 

A: The matrix H defines a eommon quadratie Lyapunov funetion for the ma- 
triees —TO and 

B: Equation (3.9) defines T = T(j:) uniquely. 

C: Control (3.8) is bounded: |u| < |. 

D: Control (3.8) brings the point ^ to 0 in time T(j:). 

E: The matrix Q is integer and even: all elements of the matrix are even inte¬ 
gers. 


Numerical experiments support the conjecture that strengthens the statement E 
of the theorem, namely, that all elements of the matrix Q are divisible by On. For 
example, for cases of one and two oscillators, the matrix has the form 


(3.10) Hi = 6 X 


1 -2 

-2 6 


1 

-9 

21 

-14 

\ 

-9 

111 

-294 

210 


21 

-294 

840 

-630 


-14 

210 

-630 

490 

/ 


4. Matching of controls 


In Section 3, we solved the problem of a local feedback control that works for 
sufficiently small initial conditions. The switching to this type of control should oc¬ 
cur at the boundary of a domain invariant with respect to the phase flow. Consider 
the invariant domains of the form 

(4.1) Ge = {y : T(r) < 6} = {y : {£2<5(0)y, (5(0)y) < 1}. 

The invariant domain G© should satisfy two conditions: 

A: The domain G© contain the “inefficiency” domain {p{x) < UC{u)} of the 
preceding control, 

B: The domain G© is contained in the strip {\Cx\ < 1 / 2 }, where G is the matrix 

(3.4). 

Condition B allows us to use at the terminal stage controls u which are less than 
1/2 in absolute value. Condition A reduces to the fact that set 1/G(A, where 
G(A, B) is a constant, is contained in G©. In other words, the following inequality 
should be fulfilled for the support functions: 

(4.2) UC{A,B)HaiD*p) < (5(0)-iq(5(0)-V,p)'/t 

where D is the matrix (3.5). It is clear that for sufficiently small U the inequality 
holds true. 

Condition B says that the value of the support function of the ellipsoid G© at 
the vector T)*“^G does not exceed 1/2 in absolute value. In other words, 

1/2 

(4.3) ((5(0)-iq<5(0)“iD*“Y,D*“Y) <1/2. 

Of course, this inequality holds for a sufficiently small 0. Once 0 is chosen, we 
have to choose the bound U for the control at the second stage in accordance with 
Inequality (4.2). Then Conditions A and B are met. The switching to the final stage 
happens upon arrival at the boundary {(n(5(0)j:, (5(0)j:) = 1 } of G©. 
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Figure 4. The structure of the suggested control for the case of 
single oscillator. 


5. Single linear osgillator 

Illustrate our preceding arguments in the simplest case of a single linear oscillator 
described by the equation (cj = 1 ) 

(5.1) 

^ ^ y = —X + R, \u\ < 1. 

We perform the separation of the phase space in three regions (see Fig. 4). 
“Basic” one is the outer region of the circle ®2 with radius 2 (circle bounded by 
solid line on Fig. 4), where the dry friction type control u = —sign(^) is applied. 
In principle, it is possible to use the of any radius r > 1 . Substantially different 
control, 

(5.2) r(x, ^) = T + 

is applied within a zone close to zero (ellipse bounded by the dot-dashed line on 
Fig. 4). Here T is the function (x^y) defined by Eq. (4.3). In our case, it has the 
form 

(5.3) - %-^24xy + %-^36x^ = 1/6. 

The close to zero zone G©, where this control is used, is the interior of the ellipse 
Q~^6y^ — Q~^24xy + 0“^36t^ = 1 with parameter 0 = 3^/^ found from the 
condition (4.3). 

The ellipse contains the disk ®a (circle bounded by the dotted line on Fig. 4) 
with radius A = (Amax)”^^^ = 0.26253..., where Amax is the largest eigenvalue of 
the quadratic form (5.3). If at the first stage we would use the circle r > 1 
instead of ® 2 , the parameter U should be A/r. 

Conglusion 

In our work a constructive method has been demonstrated for the design of a 
feedback control for a system of arbitrary large number of oscillators. In this paper, 
we have described computational aspects of the suggested approach to design of 
the control and, in particular, have shown that the design of the control in high 
energy zones reduces to solution of equation (2.3), which essentially coincides with 
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the well-studied problem of maximization a linear form on a convex hypersurface in 
As an illustrative example, the case of a single linear oscillator is considered. 
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